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We present an analysis of electron transport through two weakly coupled precision placed phos¬ 
phorus donors in silicon. In particular, we examine the (l,l)'f->'(0,2) charge transition where we 
predict a new type of current blockade driven entirely by the nuclear spin dynamics. Using this 
nuclear spin blockade mechanism we devise a protocol to readout the state of single nuclear spins 
using electron transport measurements only. We extend our model to include realistic effects such 
as Stark shifted hyperfine interactions and multi-donor clusters. In the case of multi-donor clusters 
we show how nuclear spin blockade can be alleviated allowing for low magnetic field electron spin 
measurements. 


I. INTRODUCTION 

An understanding of electron transport through multi¬ 
ple quantum dots has enabled the progression of semicon¬ 
ductor quantum information protocols from single shot 
spin readout to two-qubit logic gateJ^®!. Not only do 
transport measurements provide us with important de¬ 
tails on spin relaxation times and tunnel rates, but they 
play a vital role in aiding our understanding of the com¬ 
plex spin dynamics that occur in these systems, such as 
coherent manipulation of the electron spin^^lSl ^^d dy¬ 
namical nuclear polarisation of the Overhauser helcP. 

With recent advances in the fabrication of precision 
placed donors in silicorPtni research in this field is now 
focused on spin transport through multi-donor chains. A 
deeper understanding of the interplay between electron- 
and nuclear-spins in the dynamics of such systems is a 
prerequisite for the progression of this field. In partic¬ 
ular for the implementation of spin transport via donor 
chain^i^. So far, protocols based on spin- ch ainJI^ , CO- 
herent t unnel ing adiabatic passage (CTAPpMUIl^ gpjjj 
shuttlinJi^E^ and SWAP-gate operation^^ have been 
proposed, all of which require control of electron spin 
transport across donors. In order to further investigate 
these transport protocols it is crucial to understand the 
spin dynamics at the two donor level. However, despite 
the plethora of theoretical knowlec^ on gate-defined 
semiconductor double quantum dot^^lUM double donor 
transport has not received as much attention. 

In this paper we show how the interplay between the 
electron and nuclear spins in donor-based systems af¬ 
fects not only the charge transport, but also the spin 
transport. To understand the impact of these nuclear 
spins we use a master equation approach to conduct a 
comprehensive numerical analysis of electron transport 
through a double donor system. We investigate electron 
spin resonance (ESR) combined with Pauli spin block¬ 
ade (PSB) at low and high magnetic fields to manipulate 
and readout the electron spin states. Our most striking 
finding demonstrates that the presence of the quantised 
nuclear spin of the donors leads to a novel effect called 
nuclear spin blockade. Using this mechanism we propose 


a new spin readout protocol for the nuclear spins based 
on a measurement of the transport current. Finally, we 
analyse more realistic scenarios of inhomogeneous hyper¬ 
fine interactions across the donors, as well as the case of 
multi-donor dots, which can be shown to be immune to 
nuclear spin blockade. Throughout the paper we neglect 
the dynamical behavior of the surrounding ^®Si nuclear 
spins present in natural silicon. This interaction is much 
smaller than the donor hyperfine interactiorP^ and it has 
also been shown that Si:P devices can be fabricated in 
isotopically pure ^®Si where the absence of the ^®Si ex¬ 
tends the electron coherence time^^. 


II. TRANSPORT AT THE (1,1) TO (0,2) 
CHARGE TRANSITION 

We consider two weakly coupled phosphorus dono rs in 
a silicon lattice (approximately 15-20 nm apartpS*^ , 
and Pfl the left and right donor, respectively. Electrons 
are able to tunnel from in plane source to drain leads via 
both donors as shown schematically in Fig. Describ¬ 
ing the system is the Hamiltonian H, 

H = H,e + H,r. + Ht^+HA+Hhf, (1) 

where H^e and H^n are the electron and nuclear Zeeman 
terms, Ht^ is the tunnel coupling between the donor elec¬ 
trons, Ha is the energy detuning of the IS' 02 ) state (sin¬ 
glet state with two electrons on a single donor nuclei) and 
Hhf is the hyperhne interaction, for further details see 
Methods. Throughout the paper we refer to the Hamil¬ 
tonian in the singlet-triplet basis of the electrons with 
energies shown in Fig. By making a transformation 
from Hilbert space to Liouville space we incorporate in¬ 
coherent processes that occur during spin transportP^. In 
doing so, tunnelling that occurs from the source at the 
rate T^,, to the donors and through to the drain at the 
rate T^j, are integrated with the coherent evolution of the 
system. Using this approach we can determine the spin 
and charge dynamics of the donor system during electron 
transport. 
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FIG. 1: Electron transport through two weakly cou¬ 
pled phosphorus donors in silicon, (a) A schematic 
representation of transport through a double donor system 
during PSB. Electrons can tunnel from the source to Pl 
and from P_r to the drain at rates Pl and Ph, respectively. 
The two donor electrons are coherently tunnel coupled at 
the rate tc and have a contact hyperfine interaction with 
the nuclear spins, Al and Ar with their respective nuclei, 
(b) Eigen energies of H around A=0, between the (1,1) and 
(0,2) charge configurations at an external magnetic field of 
Bo=25 mT and with hyperfine interaction strength set to 
^i=Ai{=j4=117.53 MHz. The electron triplet states, |Tn) 
(blue), I Til) (black), and iTj^) (pink) are split by the Zee- 
man energy, 7eBo. The l^oz) (green) state is detuned from 
the ISii) (red) with an anti-crossing at A=0 due to a tunnel 
coupling set to tc = A. The (0,1) charge states are omitted 
for clarity. 


As a first demonstration of the effect of quantised nu¬ 
clear spin states we study the electron transport from 
drain to source (reverse of Fig. f^). In this scenario the 
charge cycle is (0,1)—>-(0,2)—>-(1,1)—^(O,!), where {ul, ur) 
corresponds to the electron numbers on the left and right 
donor nuclei. In the case of quantum dots it has been 
shown that the current, /qd, as a function of the de¬ 
tuning, A, is given by a known expression involving the 
coherent and incoherent tunnel rate^^, 


Iqd = 




m +{Wi^ + k) + ^^ 


( 2 ) 


where e is the electron charge and tc is the tunnel cou¬ 


pling between the two dots. However, without an ex¬ 
ternal magnetic (Bq) dependence, this equation will not 
account for any spin dynamics. In the case of donors, 
the measurable current is given by /=|e|rLP(l,l), 
where P(l,l) is the probability of being in any (1,1) 
charge configuration, including the electron triplet states 

{\Ty),\Tn)}. 

To investigate the effect the inclusion of these states 
has on electron transport we plot the difference in 
current, 51=1—Iqd, for A>0 in Fig. There is 

a peak in the current that can be seen to follow 
the point at which the electron exchange interaction, 
J=/S./2+,y (tc/2)^-|-(A/2)2, (tc=2 GHz) is equal to the 
electron Zeeman energy, 7ePo- This is analogous to the 
canonical spin-funnel experiment^. At this value of de¬ 
tuning the hyperfine interaction mixes electron singlet- 
triplet states |5 'o2)^|T'ii) via a nuclear spin flip, see 
Fig. This increase in the (1,1) state population al¬ 
lows for the electron on P/, to tunnel off to the source 
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FIG. 2: Electron transport through double donors 
from drain to source. (a) The difference in current, 
51, through a double donor as a function of detuning A, 
and external magnetic field, Bq. The yellow line maps 
J=7eBo which follows a peak in 51 due to the mixing be¬ 
tween |5o2)e->'|r{)) (contour lines guide the eye), (b) Eigen 
energies of H at Bo=250 mT as a function of the detuning, 
showing the |S'o 2 )e->'|rj)() anti-crossing. Here, Hhf allows for 
electron-nuclear spin flip-flops, increasing B(l,l). The (0,1) 
charge states are omitted for clarity. In this simulation tc=2 
GHz, Ai=AH=A=117.53 MHz, ri:,=rH=100 MHz. 
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giving rise to a larger current. In gate-defined quantum 
dots this current peak is not observed since the large nu¬ 
clear spin bath has a distribution of hyperfine strengths. 
Therefore, the position of this peak shown in Fig. 
will shift depending on the exact nuclear spin configu¬ 
ration and will be averaged out due to the fluctuating 
Overhauser fielcP^. 

For A<<0 transport across the donors becomes ener¬ 
getically unfavourable as the electrons will remain in the 
(1,1) region without first loading into the (0,2) charge 
state, see Figj^. We therefore do not expect to see mix¬ 
ing between |S'o 2 )^|T'j^) and a current peak will not be 
observed. 


III. NUCLEAR SPIN BLOCKADE 


Next, we consider the transport cycle from source to 
drain: (0,1)—>-(1,1)—>■(0,2)—^(O,!) where we can expect 
PSB due to the large energy splitting of the (0,2) electron 
singlet-triplet stated. As a consequence, if the electrons 
are in any of the (1,1) triplet states an electron on 
cannot tunnel to P^j. Therefore, tunneling can only oc¬ 
cur via the singlet states and the current is given solely 
by the |5o2) probability: /= |e| rflP(|5'o2))- 

It has been shown; however, that by applying a 
continuous-wave ESR magnetic field PSB can be lifted 
by driving transitions withi n the (1,1) electron triplet 
manifold, {|r+), |Tfi)}o|T0i)l^. This applies when the 
tunnel coupling, tc, is on the order A, hence, for the fol¬ 
lowing analysis we set tc=A. Importantly, any difference 
in the Larmor frequencies of the two electrons, dwe, will 
result in |T°^)o|S'ii) transition^. The |S'ii) state can in 
turn tunnel to liS' 02 ) via the tunnel coupling tc, allowing 
current to flow. In the presence of a homogeneous Bq 
field Su!e can only come from a difference in the nuclear 
spin orientation between the two donors. 

We examine the spin transport by initialising in the 
(0,1) charge configuration and preparing a fully mixed 
state in both the nuclear and electron spins states such 
that at t = 0, po=(|toi)(toi| + |ioi)(ioi|) ® {N} (nor¬ 
malisation omitted). We first let the system reach PSB 
after which a continuous wave ESR held is applied with 
strength Bi and on resonance at the frequency, 
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where 7 e and 7 „ are the electron and nuclear gyromag- 
netic ratios, and A]^ and Api are the nuclear hyperhne 
interactions on P^ and Pi? respectively. 

At high magnetic helds, Rq^IOO mT, there is an ini¬ 
tial peak in the current due to the tunneling of the elec¬ 
trons into the (1,1) charge region, see Fig. In a 
time ^0.1 /rs PSB is reached and the current is blocked 
(Fig. §.). After application of the ESR excitation the 
current initially spikes and coherent oscillations can be 
seen corresponding to electron spin rotations introduc¬ 
ing singlet content at a frequency determined by the 
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FIG. 3: Current during Pauli spin blockade and nu¬ 
clear spin blockade, (a) The simulated current through 
two weakly coupled donors in transport for different exter¬ 
nal magnetic fields, Bq. At t=0 the density operator is given 
by Po = (|toi)(toi| -f |4oi)(4oi|)<8) {A} (normalisation omitted). 

(b) PSB is reached in ~ 0.1 fis and the current drops to zero. 

(c) An ESR driving field is turned on after 50 ps affecting 
the spin of both donor electrons. For all values of Bo there is 
an initial spike in the current corresponding to the lifting of 
PSB followed by coherent oscillations, (d) The nuclear spin 
state projections at i3o=1000 mT are unchanged during elec¬ 
tron transport or after application of ESR due to the reduced 
electron-nuclear flip-flop events, (e) At i3o = 10 mT the nu¬ 
clear Ill'll) and lb'll) approach zero probability (<0.05 after 50 
ps), this is nuclear spin blockade. Simulations were carried 
out with tc=Ai=A_R=A=117.53 MHz, A=0, ri=r_R=100 
MHz. 
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FIG. 4: Nuclear spin readout using two weakly coupled donors in transport, (a) Schematic representation of the 
proposed scheme for single nuclear spin readout of the P/. nuclear spin showing gate pulsing, NMR and ESR driving. The 
nuclear spins are initialised in by waiting Ti. NMR can then be performed on P_l nuclei by driving at a frequency 7nBo, 
which does not effect the P_r nuclear spin due to the presence of the electron in the (0,1) charge state, (b) The charge stability 
diagram for the double donor system with colors corresponding to movements in gate space for (a). Readout is performed 
by pulsing to the (l,l)-f->'(0,2) transition and applying ESR. (c) If the nuclear spin state after NMR is then nuclear spin 
blockade does not occur, (d) However, if the spin state is then no current will flow due to nuclear spin blockade, (e) 
The current through the donors is linearly proportional to the spin down probability on Pl allowing different NMR rotations 
to be mapped using the magnitude of the current during readout. Here, i3o = 1000 mT, tc=4la=Hij=A=117.53 MHz, A=0, 
rr,=rr^=100 MHz. 


ESR field strength (7eBi ^ 140 MHz, Fig. [^). The 
current quickly finds a steady state centred around the 
oscillations from the ESR driving. During this time, the 
nuclear spins are unaffected by the electron transport. 
Figure shows the four nuclear spin state projections, 
{|M),IM),IM),IM)} = o.25. 

At low magnetic fields, i?o<100 mT, the nuclear spin 
dynamics play a larger role in determining the current 
through the donors. Here the hyperfine interactio n can 
cause electron-nuclear spin flip-flops 
This process is seen to increase the nuclear spin and 
IfllJ.) populations to ^0.3 whilst the system moves into 
PSB, see Fig. The flip-flop process is further ampli¬ 
fied when ESR is applied during which the nuclear Ifl"!).) 
and llj.'fl') populations approach zero. Consequently, this 
reduces the mixing between [S'!!) and |T®i) since 
Irrespective of the initial nuclear spin state, PSB cannot 
be lifted at low magnetic fields for double donor trans¬ 
port. We refer to this effect as “nuclear spin blockade” 
since the current through the donors is being blocked as 
a result of the nuclear spin states. 

IV. NUCLEAR SPIN READOUT 

In the same way that PSB has been used to perform 
single electron spin readout in electron transport^, we in¬ 


vestigated if nuclear spin blockade can be used to perform 
readout of a single nuclear spin. Even at high magnetic 
fields nuclear spin blockade prevents current flow if the 
nuclear spin states are aligned (lltlt) or lU-lJ-)) and allows 
current to flow if they are anti-aligned (ItlU-) or lU-lt)). 
Importantly, at high fields the nuclear spin states remain 
unaffected by the electron transport; therefore, it is pos¬ 
sible to use nuclear spin blockade as a readout mechanism 
for the nuclear spin states. 

The spin measurement protocol consists of 3 stages: 
initialisation, manipulation and readout, see Fig. |§r,b. 
We initialise in the (0,1) charge region by waiting for the 
Ti of the nuclear spins, thereby preparing them in the 
ground state, I'fl'tl'). The manipulation step involves ap¬ 
plying NMR to the P l nuclear spin (which does not have 
an electron present) at a frequency 7„i?o- Importantly, 
this will not effect the P^ nuclear spin due to the presence 
of the donor electron. The readout consists of moving to 
the (1,1)—>■(0,2) transition and applying ESR to test the 
presence of nuclear spin blockade. If the nuclear spin 
on P/, is |1J.) then current will flow through the donors 
(Fig.|§;); however, if the nuclear spin is |'f|') then the cur¬ 
rent will be blocked (Fig. W)- The current is also linearly 
proportional to the orientation of the nuclear spin; that 
is, if the total nuclear spin state is (lU-lt) -I- ItllDj/'x/S then 
the current will be half that of the |lj.'f|') state, see Fig.|^. 
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FIG. 5: Effect of inhomogeneous contact hyperfine in¬ 
teractions. At high magnetic fields where the nuclear spin 
dynamics are unaffected by the electron transport the effect 
of larger SAlr=\Al — Ab.\ reduces the current. In order 
to achieve the largest current there is an optimum Bq for 
any given SAlr (contour lines guide the eye). Here, A=0, 
rr.=rr^=100 MHz. 


V. STARK SHIFT AND DONOR CLUSTERS 



Time (|ts) Time (|ts) 


In this final section we extend the simulation work to 
examine non-idealised scenarios of Stark shifted contact 
hyperfine interactions and donor clusters (multi-donor 
quantum dots) in donor based transport, for example, 

in Ref.lini3ll36l 

The electron wavefunction around a donor nuclei can 
be distorted by electric fields which changes t he co ntact 
hyperfine interaction, known as the Stark effeclfS^^. Due 
to the stray electric fields in the silicon crystal, e.g. from 
electronic gates or charges, it is most likely that two 
donors will naturally experience different hyperfine inter¬ 
actions. This inherent Stark shift is typically a few MHz 
in frequencj®^ and will create a difference in the hyperfine 
interaction between the two donors, — Ar\. 

Figure shows the steady state current after PSB and 
ESR (as in Fig. |^) for various magnetic fields and 6Alr 
values. The current increases with the magnetic field up 
to a certain field and interestingly there is an optimum 
Bq for any given SAlr. At high magnetic fields the cur¬ 
rent decreases as SAlr becomes larger. This is due to 
SAlr decreasing Sui^. when the spins are anti-aligned, re¬ 
sulting in less mixing and a lower current. 

At low magnetic fields the effect is reversed. Here, we are 
subject to nuclear spin blockade as previously discussed, 
however increasing SAlr will only increase Soj^, resulting 
in a higher current. 

Using scanning tunneling microscopy (STM) hydrogen 
lithography it has been shown that donor clusters con¬ 
sisting of a small number of P atoms can be fabricated 
in Si that exhibit PSI^^. It is therefore possible to en¬ 
gineer a system consisting of a single donor and a two 


FIG. 6: Current through a single donor and a two- 
donor cluster, (a) Simulated electron transport for a single 
donor (Pn) and two donor cluster (2Pi{) at different mag¬ 
netic fields showing the absence of nuclear spin blockade. 
The initial density operator was po=(|toiXtoi| + lioiXXoil) ® 
{Aip/ 2 p} (normalisation omitted), (b) PSB is reached on 
a similar time scale for two single donors in transport for 
all magnetic fields, (c) After 500 rs an ESR driving field is 
applied to both electron spins. At low magnetic fields the cur¬ 
rent reaches a steady state but it does not decay to zero. The 
absence of decay in the current at low magnetic fields indicates 
that nuclear spin blockade is no longer present. At intermedi¬ 
ate magnetic fields (50 mT and 100 mT) there is a decay in the 
current; however, the final steady state current is larger than 
that obtained for two single donors in transport. In this sim¬ 
ulation, tc=A_L=A_R=117.53 MHz, A=0, Pp^Pp^lOO MHz. 


donor cluster in transport, essentially controlling the hy¬ 
perfine interaction at the atomic scale. Now, in this sys¬ 
tem the electrons will always experience a SiOg across the 
two dots, ensuring there is substantial mixing (>A) be¬ 
tween the |T°j) and [S'!!) states to lift PSB. This can be 
seen if we consider the 8 possible nuclear configurations: 

{iVip/2p}={|^)|M), l^)IW, l^)IM), l^)IM), l^)IM), 

I'D-) I til'), I'IDIM), IJDIM)}, where the first term is the 
spin state for Pp and the second term with a two nuclear 
spin state is for 2Pp. 

We simulate the current through a single donor (Pp) 
and two donor cluster (2Pp) in Fig. at different 
magnetic fields with the initial state Po=(ltoi)(toi I + 
|ioi)(4-oi|) ® {Nip/ 2 p} at t=0. With this single and two 
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donor cluster scenario, the nuclear spins take longer to 
reach a steady state configuration despite the electron 
spins reaching PSB in the same time scale as the single 
donor case, Fig. As a consequence we simulate the 
dynamics here for much longer. 

During ESR excitation the current is seen to oscillate 
with a beating due to the different hyperfine interactions 
experienced by each electron, see Fig. [§=■ For all mag¬ 
netic field strengths, the current is higher than it was in 
the single donor cases. Importantly, even at low mag¬ 
netic fields, and in contrast to the single donor dynamics 
(Fig.|§i), the current does not decay to zero. Instead, it 
finds a steady state solution centered around the coher¬ 
ent oscillations indicating that nuclear spin blockade is 
no longer present. 

VI. SUMMARY 

We have developed a numerical model to investigate 
electron spin transport through donors in silicon. Our 
findings show that there are surprising effects that arise 
due to the quantised nature of the donor nuclei. We show 
how it is possible to map out the electron exchange in¬ 
teraction using simple transport measurements and pre¬ 
dict a new current blockade mechanism as a consequence 
of nuclear spin dynamics. This is in contrast to other 
systems such as GaAs double quantum dots where large 
fluctuating nuclear spin baths averages out the quantum 
nature of the nuclear spins and these effects are not ob¬ 
served. 

Importantly, we find that for low magnetic fields, 
So <100 mT, where long electron Ti times^ and faster 
ESR rotations are possible, nuclear spin blockade pre¬ 
vents the transfer of the electron spins across the donors. 
Therefore, electron spin readout can only be performed 
at high magnetic fields; in which, the nuclear spin states 
are initialised in |f|-U.) or |U.f|-) otherwise 6oJe=0 and there 
will be no |S'ii)e-)'|T°j^) mixing. In addition, we demon¬ 
strate that there is an optimum Bq field for a given Stark 


shift across the donors to achieve the maximum current 
after ESR driving. 

Interestingly, we show how nuclear spin blockade can 
be utilised as a new readout mechanism for nuclear spins. 
This method, which only requires electron transport mea¬ 
surements and removes the need for traditional charge 
sensors will be a useful experimental tool for probing 
multi-donor interactions. The nuclear spin coherence 
times can be measured during electron shuttling w hich is 
critical for many proposed spin transport protocol^^^. 

Finally, we demonstrate that for electron spin trans¬ 
port experiments where nuclear spin blockade is undesir¬ 
able, multi-donor clusters can be used since they allow 
for lifting of PSB even at low magnetic fields. Advances 
in fabrication technologies, in particular STM lithogra¬ 
phy, offers the ability to tailor hyperfine interactions at 
the atomic scale for absolute control over the combined 
electron-nuclear spin system. 

Our results provide important insights into the com¬ 
plex spin dynamics in double donor systems. Such an 
understanding is necessary for two-qubit interaction^^ 
and spin transport protocols in multi-donor chain^Sl for 
scalable solid-state quantum computing architectureJi^. 
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Appendix A: Hamiltonian definition 

The Hamiltonian for exchange coupled donors includ¬ 
ing hyperfine, H^f and electrical detuning, Ha written 
in the singlet-triplet basis of the donor electrons is given 

by, 


Hst — Hze + Hzn + -l- Ha + Hhf 

Hze = leBo (|T+)(r+| - |Tfi)(Tfi|) + |Ho (I toi)(toi I - I ioi)(ioi I) 
Hzn = 7nBo(| I - I I) 

Fr*c = |(l*5ll)(5o2| + |^02)(^ll|) 

Ha = -A|S'o2)(5o2| 

Hnf = ArS^r’"'^ • E ■ E 

iR 


where in the Zeeman term, Ho=(0,0, Hq) and the gy- 
romagnetic ratios of the electron and nuclear spin states 


are 7e=28.024 GHz/T and 7„=17.235 MHz/T which will 






7 


define the respective resonance conditions during either 
electron or nuclear spin excitation. Throughout the ar¬ 
ticle we assume that the donors are in the weak coupling 
regime where valley-orbit effects are negligible. 

The terms and are the electron spin oper¬ 

ators for the (1,1) charge state and and lij^ are the 
nuclear spin operators for and (with donors num¬ 
bers iL=iR=i for single donors). To account for tunnel¬ 
ing out of the IS' 02 ) state to the drain we have included 
the single electron states | for) and | loi) which repre¬ 
sent the spin states of the electron in the (0,1) charge 
statd^. For these single electron states, the system has 
three particles, two nuclei and one electron spin such that 
the hyperfine coupling is only between the electron spin. 
Si, and the second donor nuclei, By considering the 
two donor nuclei quantum mechanically we can obtain 
the full system spin dynamics. 

In total our Hamiltonian is spanned by {if}(g){iV} 
where {N} represents the nuclear subspace of states 
{| 'flD')i I D'JJ'), I M), IM)} and {E} the electron subspace 
consisting of {iTn), |S'ii), IS' 02 ), | toi), I ioi)}, where i = 

+i 0, —. 

To include incoherent terms we make use of a te chnique 
commonly used in the ESR/NMR communitji^HilEIl^ 
which involves transforming from Hilbert space {Hn) to 
a higher dimensional space—Liouville space, /1„2. The 
incoherent processes can then be gathered into the dis- 

sipator superoperator, V in £„2 that can model non¬ 
trace preserving decoherence processes without the need 
of finding the operators in Hn- 


Appendix B: Rotating wave approximation 

The ESR term in the Hamiltonian (Eq. |A1[ ) is 
Hesr = COS (Ht)(s['(5 + + Si(,)), (HI) 

where we chose a realistic mT as the microwave 

magnetic field strengttP^, H its frequency and 
sj^^|.’^j,Si( 2 ;)} are the electron x-spin operators for the 
(1,1) and (0,1) charge state respectively. To remove the 
time-dependence of the ESR field we can make a transfor¬ 
mation from the laboratory frame to the rotating frame 
of the Larmor frequency (we) of the electron spins. This 
is simply a rotation about 2 : by an angle 0=LUf,t. We can 
then remove the faster rotating terms by making the ro¬ 
tating wave approximation. This transforms the Hamil¬ 
tonian according to 

Hrot=RlHRl^ -Rl^. (B2) 

This amounts to a two term addition to the Hamiltonian 

Hrot = H + Hi-F, (B3) 

where H is the non-rotating Hamiltonian, 

and F=Vla^. Hi 


contains the ESR driving terms and F is the correction 
term due to the approximation (Hi=i?®®’’/2 and cr^ is 
the 2 operator for the entire Hilbert space). 


Appendix C: Liouville Space 

To incorporate decoherence and relaxation into the 
time evolution of the density operator, a master equa¬ 
tion approach can be used. The problem is treated in a 
higher-dimensional space, known as Liouville space, £, 7 , 2 . 
The master equation in Liouville space is given bj0S, 

’^=^-t\p) + v\p) = d\p). (Cl) 

Importantly, the density operator is now a vector of 
length n^, which is operated upon by the Liouvillian 

superoperator, and the dissipator super¬ 

operator T) representing incoherent terms. When G is 
time-independent Eq. |C1| has the solution, 

|p(t))=e«5‘)|p(0)). (C2) 


Appendix D: The dissipator superoperator 


The dissipator superoperator term contains the tunnel 
rates between the source and P^, the drain and Pjj, P^ 

and Pij respectively. The elements of T) that contain 
these rates are determined by the effect they have on the 

systeirP22l. 

Two electron transport cycles were studied in this pa¬ 
per, for the hrst case: (0,1)—)-(0,2)—>(1,1)—^(O,!), the 
electrons are tunneling from the drain through the dots 
to the source. The diagonal dissipator elements in this 
scneario are. 


^[l^ii)] ~ ^il7'T+r+) 

V[\T^i)] = - Te\pto^to,) 

^[|‘5'll)] = - ^l\PSiiSii) 

B{\Tii)] = - TlIpt-t-) 

^[l-S'02)] =^i^(|ptolTol) + Iptoitoi)) 

BW toi)] =^l\Pt+t+) + + IPSiiSii)) 


h\ iol)] 


l\Ptp^tp^) + + IPSiiSii)) 

r'i? Iptoitoi) 1 


(Dl) 

where \pjk) indicates the the density operator element 
that 'D{j\ acts upon. The off-diagonal elements between 





states j and k are 


A\jk)] = ~\Pjk), (D2) 

here (r^) if j,fc=| foi) or | ioi) (I'^n), l^n)); 

otherwise, it is zero. The tunneling rates account for 
the loss of coherence between the states j and k during 
transport and must be included to ensure positivity of 
the density operator. 

The current through donors is given by the probabil¬ 
ity of the system to be in the (1,1) charge configuration 
multiplied by the tunnel rate from Pl to the drain 

I = |e|rL(|PT+T+) + \Pt°^T°J + IPt-T-) + IPSiiSii)) 

= |e|r,.(P(T+) + PiT°,) + P(rfi) + PiSn)) 

= |e|rLP(l,l). (D3) 

For the electron transport cycle: 
(0,1)—>-(1,1)—>■(0,2)—>-(0,1), the electrons are tunnel¬ 
ing from the source through the donors to the drain. So 
here the dissipator elements are essentially reversed 

^[|7’i+i)] =^|PTo.To.) 

^Il^n)] =^(IpToiToi) + IP^oaoi)) 

^[|5'ii)] =^(|PtoiToi) + IP^oaoi)) 

kT,\)]=^\pUPo^) 

2^[I-5o2)] 

h tor)] 

^[l ioi)] 

where now, for the off-diagonal elements, (Tj^) 

if j, fc=|S'o 2 ) (I toi) or I ioi)); otherwise, it is zero. 

To determine the current through the donors we cal¬ 
culate the probability for the system to be in IS' 02 ) mul¬ 
tiplied by the tunnel rate from Pfl; to the drairPSl 

I = |e|rfl|ps„.5oJ = |e|rKP(|5o2)). (D5) 

Throughout the paper we choose ri=rii;=100 MHz, 
such that it is less than A so that the hyperfine in¬ 
teraction has time to mix the electron states but still 
large enough to give an appreciable current through the 
donors. 


= - Th\ps^^So2) 

2 \PS02S02) Iptoitoi ) 

Tr 

2 \PS02S02) IP'I'OI'I'OI ) ’ 
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